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Abstract 

An uncontroversial observation of adjoint string breaking is proposed, while mea- 
suring the static potential from Wilson loops only. The overlap of the Wilson loop 
with the broken-string state is small, but non-vanishing, so that the broken-string 
groundstate can be seen if the Wilson loop is long enough. We demonstrate this 
in the context of the (2 + l)d SU (2) adjoint static potential, using an improved 
version of the Liischer-Weisz exponential variance reduction. To complete the 
picture we perform the more usual multichannel analysis with two basis states, the 
unbroken-string state and the broken-string state (two so-called gluelumps). 

As by-products, we obtain the temperature-dependent static potential measured 
from Polyakov loop correlations, and the fundamental SU(2) static potential with 
improved accuracy. Comparing the latter with the adjoint potential, we see clear 
deviations from Casimir scaling. 
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1 Motivation 



Quarks are linearly confined inside hadrons by a force called the strong interaction. 
Therefore, we cannot see single quarks: This is the basis of the stability of the matter 
we are formed of. One can study this force by analysing the energy between a static 
colour charge and a static anticharge. Unlike in the case of the electromagnetic 
force, this energy is, as a consequence of linear confinement, squeezed into a long 
flux tube. This flux tube is a string-like object. Therefore, one can ask whether this 
string actually breaks when it reaches a certain length. This breaking of the string 
corresponds to the screening of the static charges by a virtual matter-antimatter 
pair created from that very energy stored in the string. The energy of the ground- 
state of the system, the so-called static potential, completely changes its qualitative 
behaviour as a function of the distance between the two static charges and can there- 
fore be used to detect string breaking. There are two main situations where string 
breaking can be studied: (i) when one deals with static charges in the fundamental 
representation, which can only be screened by other fundamental particles, such as 
dynamical quarks or fundamental Higgs fields; (ii) when one considers static charges 
in the adjoint representation which can be screened by the gluons of the gauge field. 
To avoid the simulation of costly dynamical quarks or of Higgs fields, we simply 
consider here adjoint static charges. The bound state of a gluon and an adjoint 
static charge is called a "gluelump". Therefore, the breaking of the adjoint string 
leads to the creation of two of these gluelumps. 

Three approaches have been used to measure the static potential and study string 
breaking: 

• Correlation of Polyakov loops, at finite temperature p]. 

• Multichannel Ansatz (also known as Variational Ansatz) using two types of 
operators: one for the string-like state and one for the broken-string state 
013]. 

• Wilson loops jH E] • 

String breaking has been seen using the first two methods, but no clear signal 
has been observed using the third one. The failure of the Wilson loop method seems 
to be due to the poor overlap of the Wilson loop operator with the broken-string 
state. It has even been speculated that this overlap is exactly zero jH]. This is why 
we have a closer look at this problem, taking advantage of recent, improved tech- 
niques to measure long Wilson loops. Preliminary results have been presented in 0. 

In the next Section we recall notions about the static potential and its relation 
to the Wilson loop. In Section El we take the three methods into more detailed 
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consideration, to be able to discuss our results in all approaches. We explain the 
difficulties of measuring string breaking using Wilson loops only. In Section 01 we 
present the techniques we applied, such as adjoint smearing and improved exponen- 
tial error reduction for Wilson loops. We also discuss the methods used to analyse 
our data. For all three methods, results are shown in Section |3J followed by the 
conclusion. 

2 Static potential 

We consider a pure 577(2) gauge system with a static charge and a static anticharge 
separated by a distance R. Although we mainly focus on measuring the static poten- 
tial using Wilson loops only, we also consider the Multichannel Ansatz. Therefore, 
here we discuss the issue of measuring the groundstate and excited-state energies in 
a more general way. 

The Hilbert space of the Hamiltonian is spanned by its orthonormal eigenbasis 
ty( n \R). The corresponding energies are E < E\ < where E , the groundstate 
energy, is called the static potential. If we knew these eigenstates, we could extract 
the energies by measuring their time evolution: 

T {n \R,T) = (¥ n \R) | T T | ¥ n \R)) = e - En{R)T (1) 

where f is the transfer operator, f = J2n=o \^ n \R))e~ En ^ (^ n \R}\. Since we do 
not know these eigenstates explicitly, let us consider an arbitrary linear combination 
We can expand this state in the eigenbasis 

|$(#)) = (y {n) (R)MR)} \V {n) (R)} . (2) 

n 

The temporal evolution is given by 

| T T | $(R)) = \($(R) I ^ {n) (^))| 2 e~ E "^ T = c n e~ E ^ T . (3) 

n n 

Let us introduce a finite set of states <fri(R) which we know how to measure. We set 

$(R) = Y,*i(R)HR) ■ (4) 

i 

To be explicit, let us expand Eq.(j2J) and define the correlation matrix Vij(R,T) 

V l3 (R,T) = (4 >:j (R)\t T \4 >l (R)) . (5) 
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so that 



{$(R) | t T | = ^ aj(R)ai(R)Vij(R, T) 



(6) 



A lemma, which is proven in [Sj, states for the eigenvalues X^ n \R, T) of the 



correlation matrix V; 



\ {n \R,T) = f^(R)e~ EnT (1 + 0(e- TAEn )) 



(7) 



where f {n) {R) > and AE n = min\E n — E m \. In general, the correction term cannot 

be neglected and it will play an essential role in our study. The way to determine 
the eigenvalues \( n \R,T) and estimate the energies E n is described in Subsection 
IP1 



The finite set of states <f>i(R) is created by applying some operators on the vac- 
uum. These states are chosen to model the expected ones, the unbroken-string and 
the broken-string state. We can build a string-like state (f) s (R) by a spatial line 
S S (R) of links of length R where s denotes the number of spatial smearing itera- 
tions (see Subsection 14. or a broken-string state 4>g(R) of length R, using G(R), 
a "clover" discretisation of F^ u around the two static charges (see Subsection | 
The correlation matrix is then 



Vij(R, T) 



i(R) 



t t | MR)) 
o o\ 



hJ 



s,G 



<> 











S 8 Ss(R,T) S S G(R,T) 
GS S (R,T) GG(R,T) 



(9) 



Using the diagonalisation procedure 3 (see Subsection 14. 3|) one can reconstruct 
r (n) (i?,T) (Eq.JU) for small n, the lowest energies, hence the static potential, and 
the overlaps of the string-like state and the broken-string state with the correspond- 
ing eigenstates. In previous studies 012], this correlation matrix has been used in 
the multichannel Ansatz to show string breaking. We will confirm these results. 

Throughout this paper, we call the multichannel Ansatz the approach to obtain eigen-energies 
and -states using a multichannel approach, i.e. measuring correlations between states created by 
a finite set of operators. We call the diagonalisation procedure the numerical procedure we use 
to extract information about eigen-energies and -states from a given correlation matrix of type 
Eq.©. 
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Figure 1: The Wilson loop channel. Information about the groundstate energy, the 
so-called static potential, can be extracted. 



But since the multichannel Ansatz has been objected to (see next Section), we will 
also show that the full information about the static potential can be extracted using 
only the GG-channel or only the S S S S (R, T) channel (the Wilson loop W(R, T) (see 
Fig. |TJ)). Using the above notation, this corresponds to setting a s — 1, the other 
cii = 0. Eq.® is now given by 

(</>,(R) | f T | </>,(R)) = S S S S (R,T) = W{R,T) 

= KUR) I * (n) (^))| 2 e- En{R)T = Yl c ne~ En{R)T (10) 

n n 

We can truncate this sum of exponentials at n = I if all the states we neglect are 
strongly suppressed. For all k > I, we demand 

I ^ C fc(-ft) c -(E k (R)-E,,(m)T min (R) / n x 

ci(R) 

which defines a T m i n (R) implicitly. In particular, we consider the following two-mass 
Ansatz (1 = 1) 

W(R, T) = c e- y ^ T + Cl e~ E ^ T , T > T min (12) 

where Cq and C\ are the overlaps of the state \<ft s ) with the groundstate, which is the 
static potential, and the first excited state respectively. 



3 The three approaches 

The adjoint Wilson loop shows a good overlap with the unbroken-string but not 
with the broken-string state. This is natural: The Wilson loop observable creates 
a static quark-antiquark pair together with a flux tube joining them. Therefore, 
the broken-string state, which consists of two isolated gluelumps, has poor overlap 
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with the flux-tube state, hence with the Wilson loop. For R < Rf, the unbroken- 
string state is the groundstate, therefore in Eq.()12|) cq is large compared to c\. In 
this regime, additionally E\ > V by definition, therefore the second exponential in 
Eq.()12|) is negligible for T > T m i n . An important issue is what happens for R larger 
than Rb. The broken-string state becomes the groundstate. Its energy V(R) is 
smaller than the unbroken-string state energy Ei(R) (level crossing has occurred). 
But for small T, the unbroken-string state still dominates over the broken-string 
state because c <C c\. Of course, this domination holds only up to a temporal 
extent T = Tp, where Tp is the turning point defined by the equality of both terms 
on the right-hand side of (|12j) . For T > Tp, the broken-string groundstate becomes 
"visible" in the exponential decay of Eq. (jl2j) . The value of this Tp is crucial to be 
able to detect string breaking using Wilson loops only. 

ln^- 

Tp = E,(R) -V(R) • (13) 

Based on the heavy quark expansion, the strong coupling model of Ref. [5] leads 
to estimates of the ratio ^ and of Tp. As an example we consider the distance 
R = 12a = 1.228(1) fm, whereas string breaking occurs at R^ ~ 10a = 1.023(1) fm 
(The lattice spacing a is calculated in the beginning of Section EJ). The energy of 
the broken-string state is V(R) « 2 x M(Qg), where M(Qg) = 1.03(2)a" 1 is the 
mass of a gluelump (measured independently in Subsection 14.4(1 . Within this model, 
the ratio — is around — ~ e M (Qg) R ~ 2 x 10 5 . Therefore, the turning point is 

c o c o R=i2a 

estimated to be 

r= M(Q 9 )R 

p E 1 (R)-2M(Qg) R=i2a v ; 

using the value E±(R = 12a) = 2.34(l)a~ 1 of our results in advance. Wilson loops 
of size 12a x 42a, i.e. about 1.2fm x4.3fm, or larger are needed to observe string 
breaking according to this model. 

Based on a topological argument, Ref. jH] even suggests that there may be no 
overlap at all: c = for R > Rj,. Adding matter fields gives rise to the formation 
of holes in the world sheet of the Wilson loop, reflecting pair creation. The average 
hole size leads to two different phases of the world sheet. In the normal phase, holes 
are small and the Wilson loop still fulfills an area law, W(R,T) ~ e~ aRT , where a 
is the string tension renormalised by the small holes. This phase corresponds to the 
unbroken-string case and a screening of the static charges cannot be observed. The 
other phase is called the tearing phase, where holes of arbitrarily large size can be 
formed. As a consequence, the Wilson loop follows a perimeter law, W(R, T) ~ e~ cT . 
This corresponds to the broken-string state, since the groundstate energy remains 
constant. [H] speculates that the Wilson loop is in the normal phase, and analyticity 
prevents it from changing phase, so that string breaking cannot be seen. This would 
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explain why, for instance, the authors of [5] could not observe string breaking even 
at a distance R ~ IRb- Note however, that their temporal extent was T < 3. 

A simple argument gives a necessary condition for the observation of string break- 
ing if we only use ordinary (non-smeared) Wilson loops W(R,T). If R > Rb, where 
Rb is the string breaking distance, but T < we can relabel the R and T direc- 
tions, so that now R < Rb and T > Rb. In that case, string breaking is not visible 
since the new spatial extent is R < Rb, and the Wilson loop must still fulfill the 
area law. Therefore, both sides R and T should be larger than Rb. Replacing or- 
dinary Wilson loops by spatially smeared ones may relax this requirement somewhat. 

As already mentioned in the previous Section, an improved determination of the 
static potential can be achieved by a variational superposition of the unbroken-string 
and the broken-string states. In other words, the multichannel Ansatz enlarges the 
operator space to a multichannel approach. The unbroken-string state is realised via 
the flux tube of the Wilson loop, the broken-string state is modelled by considering 
two gluelumps - separated static charges surrounded by gluons which screen the 
"interior" colour charge. We end up with the two-channel transfer matrix of 
Eq.©. 

Nevertheless, this method has been criticised JTU]- One may claim that string 
breaking is built into the multichannel Ansatz due to the explicit inclusion of both 
states. Moreover, the behaviour in the continuum limit {(3 — > oo) must be consid- 
ered. If the off-diagonal element S S G(R,T) = GS S (R,T) is zero in this limit, string 
breaking does not actually happen: the Wilson loop does not communicate with the 
broken-string state; the eigenvalues merely cross each other at R = Rb- It is only if 
the off-diagonal elements are different from zero, that the eigenstates are a mixture 
of both states. A small overlap at different (3- values has indeed been confirmed by 
the mixing analysis in the multichannel Ansatz, as can be seen in Refs. [210], which 
deal with this issue. However, whether this overlap vanishes or not for a — > still 
has to be checked in detail. Here, we do not address this question of continuum 
extrapolation due to technical difficulties: The method we use is more efficient at 
smaller (3. Therefore we consider only one /3-value, which we choose as small as 
possible while staying in the scaling region. 

The Wilson loop method and the multichannel Ansatz work at zero tempera- 
ture, where the question to be answered is: What is the groundstate of a system 
with two static charges? In the context of the Polyakov loop method we have con- 
tributions of temperature-dependent effects, and the question to be answered is 
different: What is the free energy of a system with two static charges coupled to a 
heatbath? Nevertheless, it is an interesting issue and, as a by-product, we can also 
measure the correlations of Polyakov loops, which are in our case of adjoint charges 
-P(x) = | Tr n^o" 1 ^adj( x ' This results in a temperature-dependent potential 
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V T (R) CH: 

(P(0)P*(R)) = e - Vr{R)/r (15) 

The flattening of the potential Vr(R) at large R has been seen in QCD with 
dynamical fermions, although in practice only at temperatures close to or above the 
critical deconfmement temperature. Unlike the multichannel Ansatz, this method 
builds in no prejudices about the structure of the groundstate wave function. 

4 Technical details 

We are considering Wilson loops in the adjoint representation (for a definition see 
the following Subsection). The choice of the representation has a direct impact 
on the static potential, which is at lowest order in perturbation theory in (2+1) 
dimensions ^2] 

V P (R) ~ -cA In RA + 0{g*/A) . (16) 

C*2 is the value of the quadratic Casimir operator of the representation of the static 
charges, i.e. 

• fundamental representation: C^F) l 2x 2 = f 12x2 

• adjoint representation: C 2 (A) l 3x3 = 2 l 3x3 

The important point is that, in the regime of perturbation theory, i.e. at small 
distances R, at lowest order, the adjoint static potential V^dj(-R) differs from the 
fundamental static potential Vf UIK j(-R) by a factor |. Assuming, for simplicity, that 
the ratio remains the same at larger R (this issue is discussed in Subsection I5.5|) . the 
adjoint potential is much more difficult to measure than the fundamental one: The 
Wilson loop is W(R,T) ~ e~ v ^ T , therefore the signal decreases much faster with 
R or T in the adjoint representation. This is the price to pay if we consider adjoint 
static charges instead of fundamental static charges in order to avoid the simulation 
of dynamical quarks. We need a sophisticated method of exponential error reduction 
(see Subsection 14 .2|) to detect very small signals: The magnitude of each measured 
Wilson loop is Wi(R,T) ~ 0(1) while the average to be detected, as it will turn 
out, is W(R,T) ~ O(10 -40 ). Using ordinary methods, 10 80 measurements would be 
needed. 

4.1 Adjoint smearing 

It is very desirable to reduce contributions from excited states \I/( n ^ ) to the Wilson 
loop W(R, T): The turning point Euclidean time Tp (Eq.fJTHl)) is reduced, and the 
accuracy on the groundstate potential is greatly improved. To this end, we smear 
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adjoint links spatially. In SU(2), a matrix £/f un d in the fundamental representation 
can be mapped onto a 3 x 3 real link matrix U a d] in the adjoint representation by 

^(^fund) = 2 a ^ h ^ a % U ^M • ( 17 ) 

where the a a are the Pauli matrices; a, (3 = 1..3; i,j,k,l = 1..2. 

The smearing can be done by setting the new adjoint link as the SO(3) projection 

of the old link plus a weighted sum of the spatial staples: 

U' adi (x) = Proj 50(3) I U adi (x) + Adjoint Spatial Staple^ 
V i=i 

where we choose a = 0.5. We consider three different smearing levels: 15, 30 and 
60 iterations of Eq.ffTHj). For details and usage, see Subsection 14.2.31 We define our 
projection of an arbitrary matrix A onto A G 50 (3) by maximising TtA^A. This 
can be performed using the singular value decomposition: Every M x N matrix A 
(M > N) can be written as the product of a column-orthogonal M x N matrix U, a 
diagonal N x N matrix A with positive or zero elements (the singular values), and 
an orthogonal N x N matrix V'. 

A = UAV j (19) 

Since in our case M = N = 3, both U and V are elements of SO(3), and we get the 
projection of A onto 50(3) by 

A = Proj so(3) (I) = UV* (20) 

For completeness we also describe the fundamental smearing procedure used in 
Subsection 14.41 Instead of considering adjoint links, we use the same method but in 
the fundamental representation: 

U'{x) = Proj 5t/ ( 2 ) ^U(x) + a Spatial Staple^ (21) 

where we take a = 0.5 and project back onto SU{2) by 

U'(x) = Proj su(2) (U'(x)) = U ' {X) (22) 

; det U'(x) 



S 



4.2 Exponential error reduction 

An adjoint R by T Wilson loop in the (x, t)-plane consists of two adjoint spa- 
tial transporters of length R, which we call L(0) and "L(Ty , and two temporal 
sides of length T, which we write explicitly as U(0)t = YlJ=o '(x,y,t + i) and 

U(-R) = Ili=T-i ^idj (x + R,y,t + j) 4 . The expectation value of the Wilson loop can 
be written as 

W(R, T) = i J [DU] Tr{U(0) t L(T) t U(i?)L(0)} e~ s[u] (23) 

An exponential error reduction is possible because of the locality of the action, 
which in our case is the Wilson plaquette action. The main idea is to write the 
average of a product product of averages. 



4.2.1 Multihit-method 

One possibility to reduce the variance of the Wilson loop observable, is to reduce the 
variance of a single link contribution. The Multihit method |T3] takes the average 
of many samples of one particular link with all other links held fixed. As we will 
show now, all the temporal links 5 U®. k in Eq. (|23|) can be treated like this for 
Wilson loops with a spatial extent R > 2. In the first step we split the action as 
S = S'[U'] + J2k S'kiUk]- S'k[Uk] is the local part of the action that contains the 
fundamental link [/& corresponding to f/adj,fc- The Multihit-method can be applied 
also to Uadjj if S'{,\Uk] does not depend on U\ for k ^ I. Therefore, since we use the 
Wilson plaquette action, this condition is satisfied if R > 2. We can then apply the 
Multihit-method on all time-like links and Eq. ()23j) can be rewritten as 

[DU'} Tr (j| J dUdJUi ^' lU ^ UT^x (24) 

( f[ / dUjU^ e-w) L(0) e-^'l = 






1 / [DU] Tr ( TT Ul u ) L(T)t ( TT U adhJ I L(0) e~ s ^ (25) 



4 To make clearer the distinction between the two temporal sides of the Wilson loop, we use 
separate indices i and j for the two running time-coordinates. 

5 The situation is different considering spatial links: Since we have smeared them spatially, we 
cannot " hold all the other links hxed" . 
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where the Multihit-average is given by the one-link integral 



f dUie 



fj _ J "-^t^adj, t c 



In simple cases, as in pure SU(2), the Multihit-average can even be calculated 
analytically. Namely, S'/[Ui\ = -/3§ Tr UiW\ where W is the sum of the four (in 
3d) fundamental staples, and 

o^wjf^ , 27) 



where w = vdet W, W = W/w is the projection of the staple-sum onto SU(2) and 
J^adj represents W in the adjoint representation via Eq. ffTHj) . 

Since the variance of each time-like link entering W(R, T) is reduced, the variance 
reduction in W(R,T) is exponential in T. The coefficients have been estimated in 
[Ti] . For fundamental Wilson loops, the reduction is about (0.8 2 ) T = e _0,45T , and 
for adjoint loops about (0.5 2 ) T = e~ L39T . 



4.2.2 Multilevel-method 

Although the Multihit-method was revolutionary in 1983, the error reduction was 
not strong enough to enlarge measurable Wilson loops to temporal extents T suffi- 
cient to be able to observe string breaking. In Section El we suggested a heuristic 
argument, that T should be at least as large as the string breaking distance, which 
in our case is at Rb ~ 10a. The heavy quark expansion even results in an estimation 
of T { p st) w 42a (see Eq.flHJ)). 

M. Liischer and P. Weisz generalise the Multihit method from single time-like 
links to link-link correlators T(R,t' = na) ^3]. Using our notation from above, 

T(R, t' = na) af3lS = U*$(x, y,t + i = na) a/3 U^(x + R,y,t + j = na) lS ■ (28) 

A single Wilson loop W S (R,T) can be written, using the tensor multiplication de- 
fined by 

{T(,R, na)T(,R, (n + l)a)} afhS = T(R , TlCl) aX-ye 

T(R, (n + l)a) Xf3€5 (29) 

as 

W S (R,T) = L(0) ai {T(R,0)T(R,la)...T(R, {T - l)a)} aMS L{T)} s . (30) 

Just like in the Multihit-method where we considered the average links Ui, here 
time-slice expectation values of a link-link-correlator T(R,na), denoted by [ . ], 



10 



[T(0)][T(a)] [T(2a)][T(3a)] 



[T(2a)][T(3a)] 



;t(q)][t(o)] 



[T(3o 
[T(2a)] 
[T( 



{ 



A i A 



Figure 2: Hierarchical scheme. Using the three-level method as described in the text 
decreases the statistical error exponentially. While a one-level approach only allows 
to sample over the temporal links, a multi-level approach also makes it possible to 
update the interior spatial links, which improves the error reduction. 



can be obtained by sampling over the corresponding sublattice which is in this case 
the time-slice at time t' = na. This sublattice can be studied independently of 
the surrounding lattice provided the spatial link variables at the boundaries are 
held fixed. This is a consequence of the time-locality of the gauge action. Using a 
self-explanatory notation: 



[T(R, na)} 



1 



J sub 



[DU} suh T{R,na)e 



S[U] auh 



(31) 



the expectation value of the Wilson loop can be written in the form 



l¥(i?,T) = (L(0) a7 {[T(/?,0)][T(i?,la)]...[T(i?,(T-l)a)]}^ 75 L(T); 5 ) . (32) 

The restriction of fixed spatial links at the boundaries becomes manifest in the 
fact that only the temporal links on the time-slice at time t' = na are allowed to be 
updated when evaluating [ . ] with Monte Carlo methods. For our project, this is a 
severe obstacle which limits the efficiency of the exponential error reduction. It can 
be circumvented by using a hierarchical scheme based on identities like 



[T(R, na)][T(R, (n + l)a)] 



T(R, na)T(R, (n + l)a) 



[T(R, na)][T(R, (ra+l)a) 



(33) 



The impact on the Monte Carlo method is, that we are also allowed to sample 
over spatial links on the time-slice (n + l)a since the boundary of the so-called 

second-level average . . now consists of the spatial links in the time-slices na 

and (n + 2)a. The possibility to use this two- level scheme allows us to measure long 
Wilson loops almost up to the desired accuracy. We actually implement the following 
three-level scheme illustrated in Fig. |2l where T is restricted to be a multiple of 4: 
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W(R,T) = (L(0) Q7 {- 



[T(R, (n + 2)a)][T(R, 



[T(R, na)}[T(R, (n + l)a)} 



[n 



3)a)] 



L(r)^> 



(34) 



As a result of a rough optimisation process in the error reduction, based on 
minimising the CPU time versus the error, we choose the following parameters (see 
also The innermost averages [T(R,na)} are calculated from 10 sets of time- 

like Multihit-lmks, each obtained after rii = 10 updates. The updates alternate 
heatbath and overrelaxation steps in the proportion 1:4. The second-level averages, 

[T(R, na)} [T(R, {n + l)a)] are calculated from n<i = 160 averages of [T(R, na)} 

and [T(R, (n + l)a)], separated by 200 updates of the spatial links on time-slice 
(n + l)a. Finally, the third-level averages are calculated from n% = 165 second-level 
averages separated by 200 updates of the spatial links on time-slice (n + 2)a. 



n\ = 10 seems rather small, but can be explained using the confinement-deconfinement 
phase transition. For the ratio [T7] finds a value 



1.065(6) 



(35) 



with periodic boundary conditions. The corresponding critical temporal extent 
lattice units) for our coupling /3 = 6.0 is then 



in 



a((5 = 6.0) 



3.76 



(36) 



As a rule of thumb, one can estimate that a slice with fixed boundary conditions 
with a temporal extent of > \Nc will be "confined" (THJ. On the first level, we 
deal with time slices of extent 1, a high temperature regime, corresponding to the 
deconfined phase. Then the link-link-correlator has a large finite value, even at large 
R, and its determination is easy: n\ = 10 Multihit-averages are sufficient, and in- 
creasing ri\ further does not reduce the final error as 1/ \fn{. On the next level, time 
slices of extent 2 are in the confined phase. The signal then decreases exponentially 
at large R. We adjust ri2 so that the signal to noise ratio is about 1 for the distance 
R = 13a which we are interested in. On the third level, the situation becomes more 
complicated and the best choice of all three parameters can only be found using 
optimisation, to minimise CPU time versus error. We find that a three-level scheme 
is more efficient than a one- or two-level scheme. 
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[T(2o)][T(3a)] 



adjoint spatial transporter 




[T(0)][T(a)] 



adjoint spatial transporter 



Figure 3: The improved spatial transporters. After each calculation of second-level 
averages [T(R,na)][T(R, (n + I) a)] , we form staple-shaped transporter including 
smeared spatial links at the time-slice (n + 2) a. 



In Eq. (134)1 . we have a product of N tensors, the three-level link-link-correlators 
with a temporal extent 4a. Just like in the case of the Multihit-method, the variance 
of each tensor is reduced by a factor 8(R). Thus, the variance of the Wilson loop 
average can be reduced by as much as 5(R) N for an effort oc N5(R)~ 2 . Variance 
reduction exponential in N = T/4 is achieved. 



4.2.3 Improved spatial transporter 

Using the above technique, we are able to reduce exponentially the error coming 
from the temporal links of the Wilson loop. But there is still the intrinsic noise, 
coming from the frozen spatial links at time-slices T = mod 4, which is relatively 
large [IS]. How can we decrease it? This is what we do: To provide additional error 
reduction also in the spatial transporter, we replace the spatial transporters L(0) Q7 
and h(T)*g S with staple-shaped transporters, which are constructed in the following 
way (see Fig. E]): 

(z) After each calculation of second-level averages, denoted as I , we form 
the smeared spatial links, — , at time-slice (n + 2)a. (ii) We multiply them with 

the second-level averages to obtain the staple-shaped transporter I— I t I . (iii) 
This procedure is repeated n 3 = 165 times, each time after updating the spatial links 
on time-slice [(n+2)a] during the calculation of the third-level average. These error- 
reduced staple-shaped transporters replace the naive spatial transporters L(0) Q7 and 
h(T)% s and by contracting them with none, one, two, etc... third-level-averages one 
obtains Wilson loops at a fixed R with temporal extent T = 4,8, 12, etc... 
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4.3 Diagonalisation procedure 



Given a set of states (pi and the correlation matrix (R, T) as introduced in Sec- 
tion^ one can approximate the eigenstates correlation matrix T defined in Eq.Q, 
get information on the eigenstates \E ,( - n ) and extract the lowest energies E n using a 
diagonalisation procedure. 

For a given separation R, the correlation matrix is defined in Eq.(J5J) as 

V i] (R,T) = (<j )j (R)\f T \ ( f )t (R)}. (37) 



A naive determination of the lowest energies E n is obtained by looking for a 

A(">(fl,T) 
A(")(R,T+1) 



plateau in the ratio of eigenvalues ^RTrt^kT) °^ f° r increasing T (see 



Eq.©): 

EJR,T) = lim In - Wr ^Z ; . . (38) 

This simple method works very well, especially in the multichannel Ansatz. Never- 
theless, we want to increase the signal of the desired state as much as possible. For a 
finite basis, for small T, the eigenstates change with T. To improve T- convergence, 
we apply variational diagonalisation which consists of solving the generalised eigen- 
value problem 

V^R, T) vf\R, T, T ) = X^(R, T, T ) V^R, T ) vf\R, T, T ) , T > T (39) 

The eigenvalues \( n \R,T,T ) also fulfill Eq.© but their coefficients f {n) {R) are 
enhanced by construction, compared to the previous ones. Once the eigenvectors 
v^ n ) and eigenvalues are known, one may approximate the eigenstates as a 
superposition of the operator states using 

¥ n \R,T) = (R, T, T ) ^ {R, T, T )(j)j(R, T) = ^ a^(R, T, T )(j)j(R, T) 

3 3 

(40) 

where the constants are chosen such that the are normalised to unity. 



It is well known from the literature, that the groundstate energy can be extracted 
nicely using a single-mass Ansatz even in the broken-string regime R > Rb, starting 
with T = 0, if broken-string state operators are included in the multichannel (vari- 
ational) basis. Indeed, we can confirm this observation. The situation is different 
when one uses a pure Wilson loop operator basis, at R > Rb- Because the overlap 
of the Wilson loop with the broken-string state is so weak, one must choose To very 
large to ensure that the lowest-lying eigenstates at To and T are both broken-string 
states. Then, the high sensitivity of Eq. (JBTJj) to statistical noise renders the analysis 
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R 



Figure 4: Interacting gluelumps at distance R. Four "clovers" are stuck on a link- 
link-correlator tensor of temporal extent T. 



delicate: The matrix Vij(R,To) may not be invertible. As a trade-off, we choose a 
small value of To, To = 4a, but must use a two-mass Ansatz to account for all data 
points. 



4.4 Gluelumps 



In order to fully implement the multichannel Ansatz, we must consider the two- 
gluelump correlator, which is in the notation of Eq.® 

O <> 



<> 



GG(R,T) 



(41) 



<> 



The broken-string state can be described by the presence of two gluelumps, each 
formed by the coupling of adjoint glue to an adjoint static charge [T^j. It is of 
interest to measure the mass and the correlator of gluelumps. A simple way to 
probe the gluon field distribution, denoted as o, around the adjoint static charge 
is the so-called clover-discretisation of F^, denoted as C m>iJiU = Ci )fl „, in 
Fig.H 



I 4 

Cm,[iu = t Smeared Spatial Plaquette(m, /i, v)i 



^3-2x2 + iC^i/^ 



(42) 
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The anti-hermitian part of the clover C a approximates 

\\C m ^ v - d^J = C^ u a a = g,a 2 F, v + 0(a 4 ) (43) 

where the a a are the Pauli matrices and we average over the four oriented spatial 
plaquettes that share a corner with one end of the time-like line. We choose the 
orientation of Fig. |H (all four staples clockwise) to obtain the lowest-lying gluelump 
mass [Sj. The plaquettes are built using fundamentally smeared links (as described 
in Eq.(EI])). 

To measure the gluelump mass, one considers only one gluelump, e.g. the left 
side (Ci - C 3 ) in Fig. HJ The gauge-invariant operator CC(T) is an adjoint time- 
like line of length T using Eq. (jl7|) . A af3 (T), which is coupled at the two ends to the 
clovers (Ci and C 3 ). 

CC(T) = C«(Q)A aP (T)Cl(T) (44) 

The adjoint time-like line can be measured using the Multihit-method or by 
applying the Multilevel-idea to this particular problem. The mass M(Qg) can then 
be extracted using 

CC{T) ~ e" M(Qs)T (45) 

The overlap with the lowest-lying state is enhanced by using smeared links to build 
the clover observable Eq. (j42|) . We would like to mention at this point that the 
gluelump mass by itself has no real physical meaning, since it contains a UV- 
divergence in the continuum limit due to the self-energy of the time-like links. Only 
the difference between this divergent mass and another similarly divergent one, like 
the static potential, makes physical sense. As a consequence, the string breaking 
distance remains constant in physical units, while the energy of the level-crossing 
diverges as (3 — > oo. 

To measure the correlation of two gluelumps, one has to consider the full object 
in Fig. I3J The four clovers Ci, C4 are measured as described above. The corre- 
lation GG(R, T) of two gluelumps separated by a distance R can be measured by 
contracting the four clovers to the link-link-correlator tensors with temporal extent 
T, T(R,T). The same tensors, obtained with the Multilevel algorithm and used for 
Wilson loops, are also used here. 

GG(R, T) = C2(Q)C?(Q)T a ^ 5 (R, T)Cl{T)Ct{T) (46) 

GG(R, T) can be used in two ways: On the one hand, as mentioned in the beginning 
of this Section, it is a contribution to the multichannel matrix; on the other hand, 
we can try to extract, from it alone, the energies of the unbroken-string and of the 
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broken-string states since presumably the two-gluelump correlator has projection on 
both states. 

GG(R, T) ~ g e- v ^ T + 9l e~ E ^ T , T > T min (47) 

where V(R) is the static potential and Ei(R) the first excited state energy. The 
operator has obviously a good overlap with the broken-string state, but a poor one 
with the unbroken-string state. This situation mirrors that of the Wilson loop, de- 
scribed in the beginning of Section El For R < Rb, the unbroken-string state is 
the groundstate. Therefore go is expected to be small compared to g±, and the first 
excited state, with the energy of two gluelumps, is dominating for small T. Since 
the groundstate will be visible for large temporal extents only, we need a two-mass 
Ansatz to describe the correlator. At large distances, R > Rb, the broken-string 
state is the groundstate and also the dominating one (go gi), therefore a single- 
mass goe~ v ( R ^ T will suffice. 

In the case of Wilson loops, we attach improved spatial transporter to the link- 
link-correlators. Here, we use non-improved clovers for simplicity. Therefore, we 
have more statistical noise, which makes it difficult to extract the groundstate, if 
the turning point is large. This is the case, for distances R close to but below the 
string breaking distances Rb- For details see Subsection 15.21 

4.5 Multichannel Ansatz 

To complete the multichannel Ansatz, we must also consider the off-diagonal ele- 
ments in Eq.®, denoted S Sl G(R,T) and GS S2 (R,T). 

S S1 G(R,T) = L Sl (0) a ^T afhS (R,T)Ci(T)C s 4 (T) (48) 
GS S2 (R, T) = C'<(0)C?(0)T a ^ s (R,T)L S2 (T)} 5 (49) 

To extract their values at T = mod 4, we use the non-improved spatial trans- 
porter L S1 (0) and L S2 (T), where we smeared the links beforehand using si respec- 
tively s 2 smearing iterations. The clovers are denoted as Ci, C4. We attach them 
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to the link-link-correlators T(R,T). The complete multichannel matrix is 

/, , t ?\ 



o 



S S1 S S2 (R,T) S S1 G(R,T) 
GS S2 {R,T) GG(R,T) 



y<> <> 


<> 


V 






S15S30 




S15 G\ 


S30S15 






5*30 G 








SqqG 


\ GS15 


GS30 


GSqq 


GG J 



(50) 



This 4x4 matrix can be analysed using the diagonalisation procedure described 
in Subsection l4.31 We end up with enhanced signals for the three lowest-lying states, 
plus effective information about higher states. 



4.6 Polyakov loops 

According to Eq.([15|) we can extract a temperature-dependent potential Vj-(R). The 
correlator of two adjoint Polyakov loops can be easily built by using the link-link 
correlator tensors and the tensor-multiplication defined in Eq.([29j). 



(P{0)P*{R)) = ({■■•[ [T(R, no)} [T(R, (n + l)a) 



[T(R, (n + 2)a)][T(R, (n + 3)a)] 



} 00:77) 



(51) 



5 Results 

We are using a 3d-lattice with extent (48a) 2 x 64a at inverse coupling ft = — 6.0. 
The lattice spacing a can be obtained from the Sommer scale r 6 ^H]; which is 

6 The phenomenological interpretation of this scale is valid only for QCD. When extracting the 
lattice spacing a from the Sommer scale ro = 0.5 fm, the resulting value of a depends on the Ansatz 
chosen for the potential, and on the fitting range for the force. We chose Ansatz Ea. (|53f) . which 
includes a perturbative logarithmic term, because not including this term causes an unacceptably 
bad fit. Fitting the force over the interval 3a < R < 7a, one obtains r^ja is 4.890(1), and 
a = 0.1022(1) fm. This value changes by 0.5% under a variation of the fitting interval. Similar 
ambiguities arise when one tries to extract the lattice spacing from the string tension (era 2 ). Setting 
yfd = 440MeV w (0.45fm) _1 , one obtains a = 0.1136(1) fm, with a systematic variation of about 
0.5% with the fitting range. 
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defined by 

r^nd^o) = 1-65 . (52) 

Setting ro = 0.5 fm and comparing with the lattice result for ro/a, one obtains 
a = 0.1022(1) fm. A description of our procedure to extract the fundamental force 
-Ffund is given in Subsection 15.51 

We present our results in the following order: 

1. The static potential and excited states extracted from Wilson loops only. 

2. The static potential extracted from the two-gluelump correlator. 

3. The static potential and excited states obtained from the multichannel Ansatz. 

4. The temperature-dependent potential obtained from the Polyakov loop corre- 
lators. 

5. The comparison of fundamental and adjoint potentials and the issue of Casimir 
scaling. 

We have analysed 44 configurations, which appear to be statistically uncorrelated. 
To extract the statistical errors we apply the jackknife method. 

5.1 Wilson loops only 

We measure both the fundamental and the adjoint potential between two static 
charges. In the first case, string breaking cannot occur since the system does not 
contain particles which can screen charges in the fundamental representation. Nev- 
ertheless, we can compare our results with accurate data available in the literature 
jU |3] . We will also need these values later on, to discuss the issue of Casimir scaling 
(see Subsection 15. 5|) . In the case of adjoint static charges, string breaking should 
occur. A summary of our results is given in Fig. El where we show the fundamental 
potential multiplied by the Casimir ratio | (see Eq. ljlfij) ) and the adjoint static po- 
tential. It can clearly be seen that the adjoint static potential remains approximately 
constant for R > 10a proving string breaking at Rj, ~ 10a. The unbroken-string 
potential is also shown. 

The horizontal line at 2.06(3)a~ 1 represents twice the mass of a gluelump, whose 
evaluation is described in Subsection 14.41 This is the expected value of the static 
potential of the system, when the string is broken, since the broken-string state is 
modelled by the presence of two gluelumps whose interaction is screened. 

Excited states are not visible for the fundamental case since the shortest Wilson 
loops we consider have a minimal temporal extent of T = 4a and excited states 
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Figure 5: The adjoint and fundamental static potentials V(R) (the latter multiplied 
by the Casimir factor |) versus R using Wilson loops only. The adjoint static 
potential remains approximately constant for R> m 10a proving string breaking. 
The unbroken-string state energy is also drawn. The horizontal line at 2.06(3)a~ 1 
represents twice the mass of a gluelump. 



are already strongly suppressed. But they are clearly seen in the adjoint case for 
distances larger than the string breaking distance since the Wilson loop has very 
good overlap with the unbroken-string state which is an excited state for R > R^. 
More about excited states can be found in Subsection 15.1.21 



5.1.1 Static potential 

We start our discussion with the extraction of the fundamental static potential 
V{ nn d{R)- We consider only one level of fundamental smearing (30 iterations of 
Eq.([21)l) of the fundamental spatial links and do not consider a Wilson loop matrix 
in the sense of Eq. (J55|) . A single-mass Ansatz works nicely at all R in the temporal 
range T m ; n = 12a < T < 60a, where we have no measurable contribution of excited 
states. The extracted Vf un d(R) is in full agreement with the literature. 

From the static potential we can extract the string tension cr. This gives us a 
crosscheck with previous determinations [20] arid a way to express the lattice spacing 
in physical units. We use a string-motivated Ansatz 

V(R) ~ V + Ihx ^ + uR (53) 



20 



The Coulombic In — term follows from 3d perturbation theory (see Eg. fjlbj)) . The 
1/R term follows from the bosonic string model. 7 is a universal constant with value 
7 = |^(d — 2) in d dimensions |21j . The linear term describes confinement, and a is 
the string tension. 

We fit all parameters and find for the string tension a = 0.0625(5)aT 2 . This value 
is stable and in full agreement with (201 • Using our Ansatz Eq. (J53|) . 7 cannot be 
reliably extracted by a global fit of the static potential. Using instead 7 

' dip ' K ' 

the extracted 7 tends to the universal value 757 ~ 0.131: 7 = 0.126(12) and re- 

24 R=6a 

mains stable for R> Qa albeit with larger errors. 

In the following case of the adjoint static potential, the Ansatz Eq. (j53J) does not 
result in stable parameters, with or without the Coulombic In — term. Nevertheless, 
we include in Fig. a best fit of the adjoint unbroken-string energy in the range 
2a < R < Rb using this Ansatz. 

The extraction of the adjoint static potential works well using a single-mass 
Ansatz for T > 4 and R < Rb. At larger distances R, it is a more complicated 
matter since the string breaks and the Wilson loop has a poor overlap with the 
broken-string. This makes the two-mass Ansatz mandatory. To extract the ener- 
gies of the groundstate and first excited state, as shown in the figures, we use the 
diagonalisation procedure described in Subsection 14.31 To illustrate that the static 
potential V(R) at a fixed R> Rb cannot been determined by a single-mass, we show 
in Fig. H W(R, T) = A (0) (_R, T, T ) at R = 11a and R = 12a. Here we want to make 
use of the full Wilson loop data without distorting the ratio — of Eq. (jl2)l . We use 
a simplified version of the diagonalisation procedure to obtain X^(R, T, T ) using 
Wilson loops only: 

We have three types of staple-shaped transporter lj ,n as used in Eq. (|5T)|) . In 
the same notation, the Wilson loops correlation matrix is 

(5'i5'S'i5 Sis^o 'S'iS'S'eoX 
5 , 3o5'i5 S^oSsq S^oSqo I (55) 
SgqSis SqqSso SqoSqoJ 

1. For a fixed R, we diagonalise the matrix S Sl S S2 (R,T ), where we choose T 
so that the overlap with the desired state (e.g. the groundstate) is as large 
as possible and the signal still quite accurate. Setting T > T P is a natural 
choice. E.g. at R = 12a we choose T = 24a. 



7 We do observe an increase in 7 with increasing R, visible until R ~ 6a. This increase can be 
understood as a l/i?-correction to 7 coming from the next-to- leading term in the bosonic string 
theory [2T|- 
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Figure 6: Adjoint Wilson loop data versus T, for R — 11a and R = 12a, obtained 
from a diagonalisation procedure Eq. (J56|) . applied to Wilson loops, considering three 
different levels of smearing. A two-mass Ansatz accounts for all data points. Single- 
exponentials (dotted lines) do not. At large T, the broken-string groundstate is 
exposed. Note how small a signal can be measured. 
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2. We use the eigenvectors v (R,T ), v 1 (R,T ) and v 2 (-R, T ), where the corre- 
sponding eigenvalues fulfill X^(R,T ) > \W(R,T ) > \®(R,T ), to project 
S Sl S S2 (R,T) to the different states at all T by 

A (n) (R, T, T ) = v n , Sl (R, T ) S Sl S S2 (R, T) v n>S2 (R, T ) (56) 

\(°\R, T, To), the largest eigenvalue, contains amplified information about the 
groundstate and (R, T, T ) information about the first excited state. A^ (R, T, T ) 
is some effective value containing information about the remaining excitations. 
W(R,T) = \(°\R,T,T ) has to be analysed with a two-mass Ansatz and the ratio 

— can be extracted. 

CO 

W(R, T) = c e- v ^ T + Cl e-^ T , T > T min . (57) 

V(R) corresponds to the groundstate, E\ is the first excited state energy. The 
groundstate is the broken-string state and therefore, V(R) should be the energy of 
two gluelumps E(2Qg) « 2M(Qg). At small temporal extent T of the Wilson loop 
W(R,T) we get a larger slope than at large T, as visible in Fig. |U1 This can be 
explained as elaborated in Section |HJ At small T the signal is dominated by the 
unbroken-string state. The broken-string state can only be observed once T is large 
enough since the Wilson loop observable has a poor overlap with this groundstate. 
The ratio ^ quantifies the domination of the unbroken-string state signal versus the 
broken-string state and is related to the turning point T P (see Eq.([T3))). We obtain 
^ = 3.1 (4) 10 3 , while the prediction of the strong coupling expansion jU] was 

— m 2 x 10 5 and (see Eg . (1141) ') Tp St ^ w 42a. Our numerical determination of 

c R=12a ^ R=12a 

the turning point Tp = 22(l)a is clearly below this estimate. 

R=12a 

Note that we can detect signals down to 10~ 40 , which corresponds to 10 80 or- 
dinary measurements. Previously, only signals down to 10~ 7 have been measured, 
i.e. in a regime where the unbroken-string state is dominating over the groundstate. 
This explains why string breaking has not been observed in Wilson loops up to now. 

5.1.2 Excited states 

Excited states of the fundamental representation are suppressed too much for us to 
measure. But in the adjoint representation, we have clear information about the 
first excited state. Using the diagonalisation procedure Eq.(|55|). one source of in- 
formation is the two-mass fit of W(R, T) = \( \R, T, T ) at large distances R > R b 
where the first excited state is the unbroken-string state. Another, related, source 
of information, also for smaller R, is W\(R,T) = X^'(R,T,T ). 
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Figure 7: The static adjoint potential V(R) versus R (same as in Fig. EJ) and the 
first excited unbroken-string state energy using Wilson loops only. We also show 
the energy E^R) fEq. (joT?|) ) resulting from the relativistic Nambu string theory. The 
horizontal line at 2.06(3)a _1 represents twice the mass of a gluelump. 



In the same manner as for the static potential, we adopt here a two-mass Ansatz 

Wi(R, T) = e e- v ^ T + , T > T min . (58) 

where V(R) is the static potential and Vi(R) the energy of the first excited unbroken- 
string state (as it turns out). In Fig.0 for distances smaller than the string breaking 
distance Rf, we see a clear signal of the first excited state. However, at larger distance 
we expect contributions from at least three states: The two-gluelump groundstate, 
the unbroken-string state, and an excited unbroken- or broken-string state, since 
one expects a level-crossing of the latter two depending on the spatial distance R. 

Surprisingly, at R = 8a and R = 9a we completely miss the broken-string state. 
The explanation lies in the spectrum at R = 8a, 9a. The first three terms entering 
the expansion of the Wilson loop Eq. (jlOj) are: 

• Coe~ E ° T groundstate: lowest-lying energy state of the unbroken-string, 

• C\e~ ElT first excited state: lowest-lying energy state of the broken-string, 

• c 2 e~ E2T second excited state: first excited state of the unbroken-string. 
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Figure 8: Agreement of the static adjoint potential V(R) versus R, extracted from 
the two-gluelump correlator and the one extracted from Wilson loops only (same as 
in Fig. but shifted to the right for clarity) . The deviations at R = 8a and R = 9a 
are due to a large value of the turning point Tp, as explained in the text. 

The overlap of the Wilson loop operator with the second excited state (unbroken- 
string), C2, is larger than the overlap with the first excited state, c\. In addition, the 
difference between the two corresponding energies E\ and E 2 is small. In our case, 
the second excited state indeed dominates over the first in the accessible range of 
Euclidean times. Therefore, we miss the broken-string state at R = 8a and R = 9a. 

From our fit of the adjoint potential via Eq. (j53J) we have extracted the string 
tension a. We then consider the relativistic Nambu string [22], which predicts for 
the first excited string energy 



E^R) = \ja 2 R 2 + 2tkj(1 - — ) + const. (59) 

This relativistic bosonic string theory prediction agrees with our measured first 
excited energy remarkably well, see Fig. 

5.2 Gluelumps 

We have seen that the Wilson loop operator has an overlap with both states - 
the unbroken-string and the broken-string state. What about the two-gluelump 
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operator? At large distances R > Rb, a single- mass Ansatz is sufficient since the 
groundstate is the broken-string state, which has good overlap with the correlator of 
two gluelumps, GG(R, T). Therefore, we cannot measure any signal of the unbroken- 
string state in this regime. At distances R < Rb, GG(R, T) can be analysed using a 
two-mass Ansatz 

GG(R, T) ~ g e- v ^ T + 9l e~ E ^ T , T > T min . (60) 

At small R, the turning point Tp is small, and the two-mass Ansatz Eq. (J6Uj) works 
fine. For R = 8a and R = 9a, which approaches the string breaking distance Rb, the 
energies of the unbroken-string and broken-string state are almost degenerate, and 
the turning point value is large. As mentioned in Subsection 14 A\ we used improved 
spatial transporters to measure Wilson loops. Here, we use non-improved clovers 
and have more noise. In addition, we have only one operator-state and cannot 
apply a diagonalisation procedure. Therefore, we have difficulties to measure the 
subleading groundstate exponential decay in this regime. We lose the signal of the 
unbroken-string state before it becomes visible and cannot extract the unbroken- 
string groundstate properly. 

Fig. |S1 shows the agreement between the static potential extracted from Wilson 
loops only and that extracted from the two-gluelump correlator. This confirms that 
GG(R, T) has a non- vanishing overlap with the unbroken-string state, similar to the 
fact, that the Wilson loop has a non-vanishing overlap with the broken-string state. 

GG(R,T) can also be used in the multichannel Ansatz as explained in Section 
El It enters as a diagonal matrix element in the 4x4-matrix Vij(R,T) Eq . (|5U|) . We 
will now consider this approach. 

5.3 Multichannel Ansatz 

The full multichannel matrix Eq. (|5(J|) is obtained by including the mixing terms 
SiG(R,T) and GSj(R,T). 

5.3.1 Static potential 

We analyse the 4x4 matrix using the diagonalisation procedure Eq.()39J) at T = 4a, 
as described in Subsection 14.31 As expected, a single- mass Ansatz can be applied for 
all R using \(°\R, T, T = 4a). The results, presented in Fig. El agree with the static 
potential extracted from Wilson loops only, with improved accuracy for R > Rb. 

The three Wilson loop operator states at different smearing levels have a good 
overlap with the unbroken-string state, but a poor one with the broken-string state. 
The opposite holds for the two-gluelump operator state. To confirm this statement, 
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Figure 9: The static adjoint potentials V(R) versus R using the multichannel Ansatz. 
The agreement with the static potential extracted from Wilson loops only (same as 
in Fig. shifted to the right for clarity) is good. The unbroken-string state energy 
is also drawn (dashed line). The horizontal line at 2.06(3)a _1 represents twice the 
mass of a gluelump. 

we analyse the overlaps (R, T, T ) of Eq.(|40j). We consider the overlap of all 
three Wilson loop operator states (j = Si 5 , S 30 , S 60 ) and the overlap of the two- 
gluelump operator state (j = G) with the groundstate (n = 0) 8 . We observe for 
a ( Q ] {R,T,T ) an abrupt change from 0(1(T 3 ) (R < 10a) to 0(1) (R > 10a) and 
vice versa for (R, T, T Q ) . This indicates that string breaking actually occurs at a 
distance slightly larger than 10a. 

5.3.2 Excited states 

Considering X^(R, T, To = 4a), we get information about excited states by applying 
a two-mass Ansatz. For R > Rf,, we extract the lowest-lying unbroken-string state 
as expected. For R = 8a and R = 9a, the first excited state is the broken-string 
state. For R = 7a, the broken-string state and the excited unbroken-string state 
energies are almost degenerate. For R < 7a we extract the energy values of the 
excited unbroken-string state which are in agreement with the ones extracted using 

8 In a { " =0) (R, T, T = 4a) of Eq.® wc fix T = 8a. The results arc stable for larger T, although 
with increasing errors. 
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Figure 10: The static potential V(R), the first and the second excited states energies 
using the multichannel Ansatz. We also show the energy Ei(R) (Eq. (|59|l ) result- 
ing from the relativistic bosonic string theory. The horizontal line at 2.06(3)a~ 1 
represents twice the mass of a gluelump. 



Wilson loops only. For R < 5a we can no longer extract first excited state energies 
due to statistical noise. Finally, considering (R,T,T = 4a), for R > 7a, we 
obtain the energy of the second excited state, namely the excited unbroken-string 
state. For R < 7a we cannot extract the second excited energy due to statistical 
noise. In Fig. we show the static potential extracted from \(°\R,T,Tq = 4a), 
the first excited energy, extracted from X^'(R,T,T = 4a) and the second excited 
energy, extracted from \( 2 > (R,T,T = 4a). 



5.4 Polyakov loops 

The correlator of two adjoint Polyakov loops allows the extraction of a temperature- 
dependent potential Vr(R) (see Eq.(|15|)) 

V r (R)^-^HP(0)P*(R)) . (61) 
N t a 

The temperature of our system is T = « 30 MeV since N t = 64 is the temporal 
extent of the lattice and a = 0.1022(1) fm the lattice spacing. Since this temperature 
is quite low, contributions of excited states are negligible. Therefore, we observe a 
good match with the static potential measured using Wilson loops only (see Fig. llljl . 
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Figure 11: Polyakov loop method. The static adjoint potential V(R) versus R 
extracted from the correlator of two adjoint Polyakov loops agrees very well with 
that measured using Wilson loops only (same as in Fig. shifted to the right for 
clarity). The flattening of the potential can be observed, although we cannot extract 
the value at R = 12a due to large fluctuations. 



In the regime of string breaking R « 10a we see flattening of the potential 
indicating string breaking. For R = 12a, the signal becomes very noisy and the 
average correlator is negative. 

5.5 Casimir scaling 

Since we measure the static potentials between fundamental and between adjoint 
charges with high accuracy, we can examine the hypothesis of Casimir scaling, which 
says that the ratio of the adjoint static potential over the fundamental static poten- 
tial remains equal to the Casimir value | over a broad range of distances where both 
potentials grow more or less linearly with distance. Already in Fig. we see clear 
deviations from Casimir scaling: The fundamental static potential, rescaled by the 
Casimir factor |, agrees with the adjoint static potential at small distances R <2a 
only. This confirms earlier observations of Ref. also at /3 — 6.0, and of Ref. j2S] 
at /3 = 9.0. Therefore, we apply a more careful analysis considering forces, defined 
by 

= _V(r)-V(r-q) 
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Figure 12: Ratio of forces adJ , , as a function of the spatial separation. The 

r fund 

horizontal line at | indicates the Casimir ratio expected from perturbation theory. 
We see clear deviation at distances larger than R = 2a. The ratio seems to decrease 
linearly while increasing R. 

where rj is chosen such that the force evaluated from Eq.([62|) coincides with the 
force in the continuum at tree level [T^J. The procedure to determine the explicit 
values of rj is described in the Appendix. 

In Fig. we show the ratio ,f adJ at two different /3's. In the regime of 

^fund v**J 

perturbation theory, i.e. at small distances R, this ratio is | as expected. At larger 
distances, our (3 = 6.0 data show clear deviations. The ratio appears to decrease 
linearly with increasing distance 9 R. The less precise data at (3 = 9.0 of Ref. [23] 
seem to confirm this independence, making it unlikely to be an artifact of the lattice 
spacing. 

Our results do not necessarily contradict the work of [21], which found accurate 
Casimir scaling at large distances. Ref. [21] considers the Ad SU(3) theory, while 
we consider SU{2) in 3 dimensions. Rather, what we see might be specific to 3 
dimensions. 



9 For R > the adjoint string breaks and the force F a dj is essentially zero. Large fluctuations 
at R « Rb induce the large error on the rightmost data point in Fig. If 21 
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6 Conclusions 



We have demonstrated that string breaking can be observed, using only Wilson loops 
as observables to measure the static potential. This demonstration was performed 
in the computationally easiest setup: breaking of the adjoint string in the (2 + 1)- 
dimensional SU(2) theory. Even in this simple case, the unambiguous observation 
of string breaking, at a distance Rf, « 1 fm, required the measurement of adjoint 
Wilson loops of area in excess of 4 fm 2 , with state-of-the-art variance reduction 
techniques. A similar study in (3 + 1) dimensions with a larger gauge group will be 
challenging. 

The reason for such large loop sizes is as expected: the Wilson loop has very poor 
overlap with the broken-string. Even when the static adjoint charges are separated 
by R > Rb and the broken-string becomes the groundstate, its contribution to the 
Wilson loop area-law is subdominant. The temporal extent T of the Wilson loop 
must be increased beyond a characteristic distance, the turning point Tp, to weaken 
the unbroken-string state signal and reveal the true groundstate. We find Tp ~ 2 
fm, which explains why earlier studies, which did not use similar variance reduction 
methods, failed to detect string-breaking. While large, this turning point value stays 
well below the strong-coupling estimate of which would predict a value about 
twice as large. 

Of course, string-breaking is easy to observe, over a limited Euclidean time 
extent, if one uses a multichannel approach where a correlation matrix between 
unbroken- and broken-string states is formed and diagonalised, the latter being 
modelled by a pair of gluelumps. We reproduce in this case the results in the liter- 
ature. We also consider the two-gluelump correlator, which has poor overlap with 
the unbroken-string state, and show that the unbroken-string groundstate can be 
extracted from that correlator alone, if one allows again for a large Euclidean time 
extent. Therefore, full information about the adjoint potential is contained in each 
of the diagonal elements of the multichannel matrix. 

Finally, we looked in detail at the issue of Casimir scaling, by measuring the 

F ■ ( R) 

ratio of forces „ adJ , as a function of R. We observe clear deviations of this ratio 

•Ffund (R) 

from the perturbative value |, and an apparent linear decrease with R. A consistent 
crosscheck at a smaller lattice spacing makes this violation of Casimir scaling unlikely 
to be a lattice artifact. The situation, however, may be different in the (3 + l)d 
theory. 
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Appendix: 77 and the free scalar propagator 

We follow [12] and define forces by 

F(r ,, = _V(r)-V(r-a) m 

(X 

where rj is chosen such that the force evaluated from Eq.fjHSJ) coincides with 

F(n) = C^, (64) 

the force in the continuum at tree level. This results in 

a 

Tl = ~2n(G(r, 0) - G(r - a, 0)) ' (65) 
where G(x,y) is the massless scalar lattice propagator in 2 dimensions defined by 

-AG(x,y) — S(x,y) (66) 
We simply solve for G(x = riia, y = n 2 a) by Fourier transform, namely 

r>( \ 1 V" cos ^ Ei^i) 

G (riia, n 2 a = — 7 > — ^ , 67 

Ar2 ^ 1 EL 1 2-2cos(f/ l )' 

where N ^> ni, n 2 is the number of discretisation points, taken sufficiently large that 
G(nia, n 2 a) is known to 4-digit accuracy. Finally, we set the zero-mode contribution 
so that G(0, 0) = 0.0. A summary of our results in the range we consider is given 
in Tabled 



32 



r/a 


G(r,0) 





0.0000 


1 


-0.2500 


2 


-0.3634 


3 


-0.4303 


4 


-0.4770 


5 


-0.5129 


6 


-0.5421 


7 


-0.5668 


8 


-0.5881 


9 


-0.6069 


10 


-0.6237 


11 


-0.6389 



(r-f)/a 


n/a 


2.5 


2.3790 


3.5 


3.4080 


4.5 


4.4333 


5.5 


5.4505 


6.5 


6.4435 


7.5 


7.4721 


8.5 


8.4657 


9.5 


9.4735 



Table 1: (left) The scalar lattice propagator G(r,0) in 2d. (right) The naive deriva- 
tion points (r — |) compared with 77 in the range we consider. 
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